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A class of topological superconductors in two dimensions support gapless Majorana edge modes 
that are protected by time-reversal symmetry. Here we study the quantum phase transition (QPT) 
associated with the destruction of these edge modes as magnetic order develops spontaneously along 
the edge. It was recently conjectured that this quantum critical point has space-time supersymme- 
try, that emerges dynamically by tuning only one parameter. Here, we put this assertion to test by 
constructing a one-dimensional lattice model of the edge, in which the role of time reversal is played 
by a non-local symmetry. We solve this model using state-of-the-art Density Matrix Renormaliza- 
tion Group (DMRG) method, and indeed find emergent supersymmetry at the symmetry breaking 
transition. In particular, we find that the QPT lies in the tricritical Ising universality, even though 
it is generically accessible by tuning only one parameter. We provide a detailed description of the 
QPT, and discuss the consequences of the emergent supersymmetry. 

PACS numbers: 



Contents 



I. Introduction [T] 

II. Magnetic QPT at the boundary of 

topological superconductor [2] 

A. Field theory of the transition [5] 

B. A lattice model for the transition [5] 

III. Numerical solution via DMRG and the 
emergent supersymmetry S] 

IV. Conclusions [S] 
References [6] 



I. INTRODUCTION 

In this paper, we solve the problem of magnetic phase 
transition at the edge of a two dimensional time-reversal 
invariant topological superconductor (TSC)^ ^. We ac- 
complish this task by constructing a lattice model that 
captures the physics of this quantum phase transition 
(QPT)[3 and we solve the model using state-of-the-art 
Density Matrix Renormalization Group (DMRG)^. We 
find that the QPT has emergent space-time supersymme- 
try and lies in the tricritical Ising universality clas^^^^^, 
even though it is accessed by tuning only one parame- 
ter. This is in complete agreement with the recent field- 
theoretic arguments that addressed the same problenP. 

The discovery of Z2 topological insulators'^"^ has led 
to an explosion of activity in the field of topological 
phases. A complete classification of non-interacting topo- 
logical phases ^^ f^^ l revea ls a class of superconductors 
(class Dili TSC)P2lIlIl9l, which supports Maj orana ed ge 
states protected by time-reversal symmetrj^E^I^nHIll 
The well known B phase of superfiuid He-S^S is a re- 



alization of a topological phase in this class , and possi- 
ble solid state realizations have also been proposed^^"^. 
The Majorana edge modes of Dili TSC are the supercon- 
ducting analogs of the celebrated Dirac fermion modes in 
topological insulators (TIsjP^. Since these modes are 
symmetry-protected, spontaneous breaking of the time- 
reversal symmetry provides a natural route to gap them 
out, without closing the fermionic energy gap in the bulk. 
For example, electron correlations often lead to magnetic 
order, which breaks time reversal symmetry. A natural 
question is: how do the edge modes evolve as the mag- 
netic order sets in? 

In a recent woriP, it was argued that the QPT where 
the magnet order develops along the boundary of a 
Dili TSC has emergent space-time supersymmetry. The 
field theory for the QPT consisted of gapless bosonic 
modes of the critical magnetic fiuctuations strongly cou- 
pled with the boundary Majorana edge states, such that 
the bosonic modes are superpartners of the Majorana 
fermions. Specifically, in d=2, it was argued that this 
is the tricritical Ising transitioiJ^^^ accessed by tuning a 
single parameter. In this paper, we construct an explicit 
one-dimensional lattice model which, in the continuum 
limit, reproduces the theory of the boundary of a TSC 
coupled to the magnetic order parameter. We simulate 
this lattice model using state-of-the art density matrix 
renormalization group (DMRG) and find unambiguous 
verification of the claim in Ref.l^l that the QPT corre- 
sponds to (supersymmetric) Ising tricritical point and is 
accessed by tuning only one parameter. We note that 
other constructions of supersymmetric condensed mat- 
ter models have appeared in the li teratu re, but they all 
require tuning multiple parameter^^Hini -^^e emphasize 
that our situation is different, in that it exploits the spe- 
cial properties of the boundary states of topological su- 
perconductors, to access a supersymmetric state on tun- 
ing a single parameter. 

A crucial property of our one-dimensional lattice model 
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is that, in the magnetically disordered phase, it supports 
gapless helical Majorana edge states of two-dimensional 
TSC. This in itself is an interesting feature because one 
might naively expect that the helical Majorana modes 
cannot be obtained in a strictly one-dimensional model 
without fine tuning. We achieve this task by observing 
that the free helical Majorana modes in one-dimension 
are d ual to the critical point of a transverse field Ising 
m odeP^ *^ which has an additional self-duality symme- 
ir?/2lMl Our model retains this duality symmetry even 
after coupling to the magnetic degrees of freedom, which 
leads to protected helical Majorana modes without fine 
tuning, as long as the magnetic moments are in the dis- 
ordered phase. 
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Figure 1: The unit cell of our lattice model consists of two 
Majorana fermions (x and x) s-^d two spin-| spins (na and 
fib) or equivalently, three spin-| spins (a,iJ,a and fj,f,). The in- 
teractions between the spins are specified by the Hamiltonian 
H in the text (Eqn|5|. 



II. MAGNETIC QPT AT THE BOUNDARY OF 
TOPOLOGICAL SUPERCONDUCTOR 

A. Field theory of the transition 

In this section we introduce the problem of the mag- 
netic instability at the edge of a two dimensional DIII- 
TSCP. Owing to the topological band structure of a TSC, 
there exist "helical Majorana modes" at its boundarjW2l. 
This means that the edge states consist of non-interacting 
Majorana fermions XR^Xl that move in the opposite di- 
rections and are described by the following quantum ac- 
tion: 



•Sx y" drdx [xB. {dr + idx) XR. + XL {dr - idx) Xl] 



(1) 



where we have set the velocity of the Majorana 
fermions to unity. The above action has an inbuilt 
time-reversal symmetry T which acts on the Majorana 
fermions in the following way: 
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. In quantum mechanics, the time-reversal operation is 
anti-unitary that is, it acts on the imaginary number 

= i as i —i. This implies that a potential mass 
term Hgap = J dx ixlXr is not allowed in the action 
since it changes sign under T. This is the reason that 
the Majorana modes at the edge of a TSC stay gapless 
as long as the time-reversal symmetry is not broken. 

We are interested in studying potential instabilities of 
the edge that gap out the Majorana modes. The mag- 
netic fluctuations can be described via a conventional 
theory for an Ising variable (j): 



drdx 



(2) 



Under time-reversal (f) reverses its sign and therefore the 
coupling between the Majorana modes and the magnetic 
fluctuations, to the leading order, is given by: 



Sint = / drdx xlXr < 



(3) 



Sint is clearly invariant under T. The full action is given 

by 



S — Sy S(h 



Si', 



(4) 



Clearly, when Ising fleld </> orders, it spontaneously gener- 
ates a mass gap cx g{(j)) for the Majorana fermions. Our 
aim is to determine the critical properties of the corre- 
sponding QPT. 



B. A lattice model for the transition 

We now construct a novel one-dimensional lattice 
model that reproduces the physics of the magnetic phase 
transition described in the previous subsection. Naively, 
one might think that the edge of a two dimensional topo- 
logical phase may not be realizable faithfully by a strictly 
one-dimensional model, and therefore one may need to 
simulate a full two-dimensional model to describe the 
QPT. We achieve this seemingly difficult task by ex- 
ploiting the self-duality synimetry of the critical 1+1-d 
transverse fleld Ising modePSES which is also dual to 
the gapless helical Majorana modes via Jordan- Wigner 
transformation^^l^Sl xhe self-duality symmetry simu- 
lates the aforementioned time-reversal symmetry T in 
a non-local manner. 

The unit cell of our lattice model (FigjT|) consists of two 
Majorana fermions and two spin-^ spins: Xi^XiJ fli,a, fii,b 
or equivalently, three spin-^ spins : <Ti, jli^a, f^i.b where 
the labels a and b serve to distinguish the two different 
fl spins within the unit cell. Similar to the action S in 
Eqn^ the Hamiltonian of our model consists of three 
terms: 



H = Hi + H2 + H3 



(5) 
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where 



i 

i 

H2 = - Yl i^^lafJ■lh + ^J■lb^^^+l,a + Hpla + ^^lb)) (8) 

i 

H3 ^ K"(mL + - + Ai'+l,a)] (9) 

i 

= {X-iXiifJ-la + l^lb) - X^X^+l{^J■lb + M?+l,a)) 

i 

Above, the Eqns|7] and [TO] follow from Eqnsj6] and [9] 
via the well-k nown Jordan Wigner transformation on the 
spin variableJ^^EH. 

One can also write down a slightly simpler model where 
is replaced by 

^3 = ig^ixiXti^la-XtXi+iiJ-lb) (10) 

i 

i 

Above, Hi is the transverse field Ising model at criti- 
cality for the a spins which is dual to the free Majorana 
fermions (Eqnj7| via Jordan- Wigner transformatioiP^*^. 
Therefore, its continuum limit corresponds to the helical 
Majorana mode action in EqnlT] H2 corresponds to 
the transverse field Ising model for the /i spins at a trans- 
verse field h. In the continuum limit, the corresponding 
action is given by (Eqnj2]). Finally, H3 describes the 
coupling between the two Ising models and precisely re- 
produces the term Sint in the continuum limit. The pa- 
rameter h is the analog of the r term in action S in Eqn|4] 
while we retain the notation g for the coupling between 
the two Ising models as before. 

One might naively expect that in our model, the a 
spins would generically be off-critical since after coupling 
to the /2 spins, the coefficient of — erf cr|^_j^ term in Hi 
may renormalize and become different than that of the 
term —J^i'^f- This will spoil the correspondence with 
the helical Majorana modes which stay gapless as long 
as scalar field (p does not condense. Interestingly, the 
above models (with either H3 or H!^ term) have an intri- 
cate symmetry, which we denote by T', that ensures that 
a spins stay gapless throughout the whole phase when fl 
spins are disordered. In particular, the above model is in- 
variant under: afaf^^ A <7f_^^, af ^ ^^la ^ 

Physically, this symmetry exploits the fact that the 
transverse field Ising model is self-dual at the criticality. 
The self-duality interchanges the terms — J^i '^t'^t+i ^^d 
— ^ • erf in Hi and this is precisely the action of the sym- 
metry T' on the (T spins. Since the symmetry is respected 
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Figure 2: The evolution of the central charge for g — 0.5 
as the parameter h is varied in Hamiltonian H (Eqn[5|. The 
central charge jumps from c = 0toc = 7/10as one approaches 
the transition between the two phases (the phase diagram is 
shown in Fig|4|. It settles down to c = 1/2 on the other side of 
the transition since the corresponding phase is described by 
the critical point of the 1+1-d transverse field Ising model. 
We note that the jump from c = to c = 7/10 becomes 
sharper as the total system size increases, since the finite size 
effects decrease. 



by the full interacting Hamiltonian H = H1+H2+H3 (or 
H' = Hi + H2 + H'^), the spins a remain gapless unless 
the T' symmetry is broken spontaneously. This happens 
when (fiz) 7^ 0, in the exact analogy with the action S in 
Eqn|4] where (</>) provides mass gap to the helical Majo- 
rana modes. Thus, for a given coupling 5, as one varies h 
from large to small values, one begins in the phase where 
the spins are disordered, implying a gapless phase with 
c — 1/2. On lowering h, there is an ordering transition 
at hc{g), below which the fi spins order, i.e. (/i^) 7^ 0. 
This leads to the gapped, symmetry broken phase. The 
critical point at hc{g), as in Figure |4] is the primary fo- 
cus of this work. We note that the unusual symmetry of 
the model under study is refiected in the fact that the 
disordered state of the spin model has a three fold de- 
generate ground state, since one of the two degenerate 
ordered states of corresponds to the ordered phase of 
the cr^, which is itself two fold degeneratd^. Thus, in 
this respect it resembles the first order line between the 
disordered and and ordered phases in Ref^. 

Having argued that our model reproduces the physics 
of the action S (Eqn|4|, universality of the critical phe- 
nomena implies that it can also be used to study the 
critical properties of the magnetic phase transition, to 
which we now turn. 
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III. NUMERICAL SOLUTION VIA DMRG AND 
THE EMERGENT SUPERSYMMETRY 
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Figure 3: The scaling of the von-Neumann entanglement en- 
tropy and the Renyi entropies ^2 and S3 at the critical point 
between the two phases in Fig|4] The scaling of the entangle- 
ment entropy fits well to the predictions from analytical cal- 
culations (Eqn |12[ ) with a central charge c — 7/10, confirming 
that the transition lies in the tricritical Ising universality. 
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Figure 4: The phase diagram of the model in EqnjS] obtained 
using the DMRG method. The Ising critical phase is identified 
by a gapless CFT with central charge c = 1/2 and disorderd 
fl spins, while the gapped phase corresponds to ordering of 
fl spins which interact with the a spins via the term in 
Eqn[5]to gap them out. The line separating the two phases 
correspond to tricritical Ising CFT with central charge c = 
7/10. The central charge along the whole g = line for h ^ 1 
is 1/2 since the a and fl spins decouple along this line, while 
the point 3 = 0, /i = 1 corresponds to a multicritical point 
with central charge unity since the fi spins are also at the 
Ising critical point. 



Wc solve the Hamiltonian H in Eqnj5] numerically us- 
ing DMRG"* calculations. We simulate lattice sizes upto 
140 unit cells with periodic boundary conditioners] and 
keep up to 1024 states to ensure the convergence of 
the entanglement entropy in DMRG optimization. For 
smaller system with 8 unit cells, we have checked that 
our DMRG results are identical to exact diagonalization 
results. The DMRG is a very useful tool for our purposes 
since it provides a direct access to the central charge of a 
theory, which uniquely determines the critical theory in 
our model. To this end, we note that one of the primary 
outputs from a DMRG calculation is the von-Neumann 
entanglement entropy S{1) for a bipartition of the sys- 
tem of size L into two subsystems of sizes / and L — I. 
For a subsystem A, it is defined as S{A) — — Tr pA log pA 
where pA is the reduced density matrix for the subsystem 
A obtained by tracing out the rest of the system. For a 
conformal field theory (CFT), S{1) is known to scale ad^ 



0(1) 



(12) 



where c is the central charge of the CFT and 0(1) 
denotes subleading corrections. One observes that the 
above expression is symmetric under I L — I ai the 
leading order which is a consequence of the fact that 
S{A) = S{A) for an arbitrary subregion A. For a gapped 
phase, c = and the entanglement entropy saturates to 
a constant value as I increases. One may also look at 
higher moments of the reduced density matrix, which go 
by the name of "Renyi entropies" 5„ = logTr p^. 

For a CFT, S'„(0 scales similar to S{1) in Eqniljp, 



i)log(^sin(^ 



0(1) (13) 



One may easily show that all the central charges in- 
volved in the problem of our interest (EqnjSj) are less than 
unity ^. and this uniquely determines the corresponding 

Figj2] shows the evolution of the central charge for 
g — 0.5 as the parameter h is varied. One observes that 
the central charge stays at zero at small h indicating a 
gapped phase. It jumps to « 0.7 at he ~ 1.62 and then 
settles down to « 0.5 at ft, > 1.62 side. We further note 
that the jump from c = to c = 7/10 becomes sharper 
as the system size increases from L = 100 to 140 with 
the decrease of the finite size effect. Fig jS] shows the en- 
tanglement entropy as a function of the subsystem size 
I at the critical point between the two phases, that is, 
'At g — 0.5 and he — 1.62. We find that von-Neumann 



entropy iS'(^) fits perfectly to Eqn 12 to yield a central 
charge c = ^ at the critical point. Furthermore using 
c = by fitting S{1), we demonstrate the agreement 
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via Yukawa interactiorP^l (Eqnj3| (the reader may note 
that in supersymmetric notation, the doublet of fermions 
X — [X-R Xl]^ counts as a single fermion). Furthermore, 
for massless fermions and bosons, the = 1 supersym- 
metric Wess-Zumino model is known to be equivalent to 
the tricritical Ising model. To see how supersymmetry 
may possibly be realized by the action S in Eqn|4] let 
us set the bare value of the couplings r = 0,c — 1 and 

2 

further assume that the relation u = ^ holds . If so, one 
may easily verify that S is invariant under the following 
'fermionic rotation!^ 



Sx 



ex 



(14) 
(15) 



Figure 5: The top figure shows the scahng of the correiation 
function (/i^/i^/) while the bottom one shows the bond-bond 
correlation function (/i^/ir+iMr'/^r'+i) ~ (MrA'r+i)^ (^^ denote 
the bond operator B^. = in the label of the figure) for 

the Hamiltonian H (Eqnj5| at the critical point between the 
two phases in Fig|4] Emergent supersymmetry at the critical 
point implies that the difference in the power-law exponents 
for these two correlation functions is exactly two while the 
precise values of the exponents themselves are also found to 
be consistent with the tricritical Ising model (see the main 
text for details). 



Here e = [ei £2]"^ is an arbitrary two-component Majo- 
rana variable, i = eiay, x = [XR Xl], It = cry, = cTx, 
and the Pauli matrices a act on the spinor index of 
fermions x ^^'^ The numerical solution of our lattice 
model demonstrates that J\f — 1 supersymmetry emerges 

dynamically at the QPT by tuning only one parameter 

2 

(in contrast to the fine tuning of four parameters: r, c, ^ 
and the Majorana fermion mass). Supersymmetry im- 
plies that the velocity as well as the anomalous dimen- 
sion of the Majorana fermion x equals that of the scalar 
(/>, a direct consequence of invariance under 'fermionic 



rotation' , Eqns 14 and 15 In particular 



between numerical results and CFT predicted Renyi en- 
tropy Sn in the inset of the Fig|3] The central charge 
c = jq between the Ising critical phase and the gapped 
phase implies that the critical point associated with the 
magnetic ordering corresponds to tricritical Ising theory, 
in agreement with the field theoretic arguments in Ref.l^. 

By repeating the same procedure as above for different 
values of the parameter g, one can obtain the whole phase 
diagram in the g — h plane (FigQ . Consistent with our 
prior assertions, our lattice model exhibits a whole phase 
at larger h side that is Ising critical without any fine- 
tuning and thus faithfully captures the edge theory of a 
two-dimensional TSC. This phase is separated by a line 
in the g — h plane that describes the phase transition 
between the Ising critical phase and a fully gapped phase 
that breaks the T' symmetry spontaneously due to the 
condensation of /i^ spins. Most importantly, we verify 
the scaling of entanglement entropy, Eqn |12| with c = 
j^, along the whole critical line that separates the two 
phases. This provides an unambiguous evidence that the 
critical theory is given by the tricritical Ising model. 
Emergent Supersymmetry and its Consequences: 
The tricritical Ising model in two dimensions is rather 
special: it has space-time supersymmetry"^-^^^. Indeed, 
the form of the action S in Eqn|4] is ready-made for 
A/" = 1 supersymmetry in 1+1-d, which is the theory 
of a single real fermion coupled to a single real scalar 



(XL(r)xL(O)) ^ (Xi?(OXi?(0)) ^ ^ (16) 

(</'M'^(o)> - (17) 

with 771 = ?72 = 0.4. Supersymmetry also implies that the 
difference of the scaling dimension of the operators /if 
and fiffif^^ is precisely unitjEH. SpecificalljESl, (/if Mr') ~ 

\r-r'\0-^ while (/if/if+i/if,/if, + i) - (/if/if+i)2 - |^_^\|2.4 . 

We verified this particular prediction in our DMRG sim- 
ulations and the results are shown in Fig [5] 

IV. CONCLUSIONS 

In this paper, we studied the quantum phase transi- 
tion (QPT) associated with the magnetic instability of 
the boundary states in two dimensional Dili class topo- 
logical superconductors (TSC) . We found that this single- 
parameter tuned transition has emergent space-time su- 
persymmetry at the critical point and lies in the tricrit- 
ical Ising universality class, in agreement with the re- 
cent field-theoretic arguments. To derive this assertion, 
we constructed a concrete lattice model for the phase 
transition which we solved using the DMRG method. A 
physical scenario leading to such a transition can arise 
in several ways. One can imagine a situation where the 
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material has two bands, one of them contributing itin- 
erant electrons that form a topological superconductor 
while the other band consists of localized moments that 
provide the Ising magnetic variable. Alternatively, one 
can also dope the system with magnetic atoms in a uni- 
form manner so that the disorder effects can be ignored. 
In the latter scenario, preferentially doping the surface 
could trigger a surface transition while the bulk retains 
time reversal symmetry. 

We note that from a field theoretic point of view 
a supersymmetric fixed point is generically multicriti- 
cal. For example, the action S is supersymmetric when 

2 

r = 0,M=^,c=l and the Majorana fermion mass is 
also set to zero. Indeed, there already exist lattice mod- 
els in literature that support supersymmetric multicriti- 
cal point^2IlS What is special about our model is that 
it leads to supersymmetry by tuning only one parameter 
due to its intrinsic dynamics. We note that higher di- 
mension al an alogs of similar phenomena have also been 
proposecPI^^lBZl and it would be interesting to realize lat- 
tice model in higher dimensions that can be numerically 
simulated and exhibit the phenomena of emergent super- 



symmetry. 

The emergent supersymmetry at the criticality has 
many remarkable consequences. We already mentioned 
that it leads to the equality of the anomalous scaling di- 
mension of the Majorana fermion (which can be extracted 
using ARPES experiments) and the anomalous dimen- 
sion of the magnetic order parameter (measured using 
neutron scattering experiments). It also implies emer- 
gent Lorentz invariance at the criticality, which leads to 
the equality of the speed of the Majorana fermion and the 
critical fluctuations associated with the magnetic order 
parameter. Finally, we note that one can also imagine 
situations where the magnetic order sets throughout the 
bulk of the TSC; and we refer the reader to a separate 
analysis of such transition^. 
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We thank David Huse for discussion on this point. 
Strictly, the duality from fermions to the Ising chain dic- 
tates specific boundary conditions, which however are 
unimportant for the physical quantities we study here. 
Hence we always work with periodic boundary conditions. 



